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Abstract. Quantum mechanics must be regarded as open systems. On one hand, this is due to 
the fact that, like in classical physics, any realistic system is subjected to a coupling to an uncon- 
trollable environment which influences it in a non-negligible way. The theory of open quantum 
systems thus play a major role in many applications of quantum physics since perfect isolation 
of quantum system is not possible and since a complete microscopic description or control of the 
environment degrees of freedom is not feasible or only partially so [1] . Practical considerations 
therefore force one to seek for a simpler, effectively probabilistic description in terms of an open 
system. There is a close physical and mathematical connection between the evolution of an 
open system, the state changes induced by quantum measurements, and the classical notion of a 
stochastic process. The paper provides a bibliographic review of this interrelations, it shows the 
mathematical equivalence between markovian master equation and generalized piecewise deter- 
ministic processes [1] and it introduces the open system in an open observed environment model. 



1. Introduction 

The dynamics of open quantum systems plays a central role in a wide class of 
physical systems. Usually, the dynamics of an open system is described in terms 
of the reduced density matrix p$ (t) which is defined by the trace over the environ- 
ment degrees of freedom [1, 2, 3]. On the ground of the weak-coupling assumption 
and the Rotating Wave Approximation the dynamics may be formulated in terms 
of a quantum dynamical semigroups which yields a Markovian Master Equation 
[4]. However, the dynamical equation, thus obtained, is very often untractable. 
This fact has encouraged some physicists [1, 5, 6, 7] to look for alternative ways 
to describe open systems. Instead of representing the dynamics of an open system 
by a quantum Master Equation for its density matrix, it is formulated in terms 
of a stochastic process for the open system's wave function. The stochastic repre- 
sentation of quantum Markov processes already appeared in a fundamental paper 
by Davies [8] and was applied [9] to derive the photocounting formula. While the 
theory was originally formulated in terms of a stochastic process for the reduced 
density matrix, in the last decade it has been proposed [1, 5, 6, 7] as a stochastic 

I thank Prof.N. Messina and Dr.A.Napoli to introduce me on the subject and to give me the 
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evolution of the state vector in the reduced Hilbert space (for a review see [10]). At 
the same time Carmichael has developed the idea of an unravelling of the master 
equation in terms of an ensemble of quantum trajectories. His theory is applicable 
only to a particular class of quantum systems (the photoemissive sources) and it in- 
duces to think that this treatment is equivalent to the Master Equation approach. 
In the last few years F.Petruccione and H.P.Breuer [1] generalize and provide a 
mathematical formulation of Carmichael's idea of quantum trajectories [5] and of 
the Monte Carlo wave function method [6]. In this contest their main result is 
to demonstrate that the dynamics given by the most general Master Equation in 
Lindblad form can be represented as a piecewise deterministic process (PDP) [1] 
ip{t) in the Hilbert space of the open system. The physical basis to achieve their 
aim is provided by continuous measurement theory [1]. 

The link [11] between the first and second way of describing open quantum 
systems is, however, only one way in the sense I am going to explain. The cited 
paper contains a flow diagram that is a clear and straightforward picture to review 
Concept and methods in the theory of open quantum systems, as the same authors 
title their work. The diagram shows us that PDPs, describing a selective level 
of measurement, implies the non- selective one, appearing in the form of a Lind- 
blad markovian master equation, while the opposite has not been demonstrated. 
My aim is to put an arrow in the opposite direction in order to demonstrate the 
equivalence between the two approaches to open systems dynamics, at least under 
Born-Markov approximation. To achieve my goal I start from microscopic mod- 
els and, exploiting the same approximation leading to the most general Master 
Equation (not in Lindblad form 1 ), I solve this last exactly (at a bath temperature 
T = 0) obtaining an operatorial expression for ps(t). Moreover, in the contest 
of optical quantum system I derive an expression for ps(t) that is in accordance 
with the Carmichael's one, but, differently from that, the mine is applicable also 
when the Master equation is not in the Lindblad form in which cases Carmichael's 
solution is very often merely formal, as the same author underline in recent papers 
[12, 13]. My mathematical tool, here called NuD Theorem, is applicable to a wide 
class of systems, provided that they satisfy the hypothesis necessary to make it 
working. Optical physical systems well satisfy the conditions of NuD theorem's va- 
lidity and in this contest I have studied, as example of monopartite system, a single 
mode cavity, and as example of multipartite systems, two two-level dipole-dipole 
interacting atoms and iV two-level not-directly-interacting atoms placed in fixed 
arbitrary point inside a loss cavity [15, 16]. In the investigation of the monopar- 

1 Evcry Markovian Master Equation can be put in Lindblad form and this, in general, introduce 
simplification in the further calculations, but because of the difficulty to recast the equation in 
this form the results obtained are in general merely formal. For this reason I prefer to work with 
non-diagonal Markovian Master Equation, by virtue of which my solution is applicable to the 
more of the known open-system theoretical models -from matter-radiation interaction model to 
the spin-boson model-, more easily than other approaches. 
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tite systems the novelty is constituted by the operatorial way to approach to the 
master equations of the systems, being already known in literature [1, 2, 3, 17] 
their dynamical properties. I reproduce for example the photocounting formula in 
order to appreciate the easyness of application of my method. Moreover I choose 
to analyze monopartite systems at a temperature T = to highlight that, differ- 
ently to multipartite ones, spontaneous emission provokes decoherence phenomena 
that, inevitably, guides the first kind of systems to their ground states. On the 
contrary, multipartite systems can exhibit collective properties induced by the 
common reservoir [15, 16, 18, 19, 20, 21, 22, 23]. This general feature already 
appeared in some fundamental papers [24, 18] in which the interaction between 
atomic dipoles, induced by electromagnetic field, could cause the decay of the 
multiatom system with two significantly different spontaneous emission rates, one 
enhanced and the other reduced. 

In two recent papers, as an example of multipartite system, I have investigated 
the dynamics of a couple of spontaneously emitting two-level atoms, taking into 
account from the very beginning their dipole-dipole interaction and and N two- 
level not-directly-interacting atoms placed in fixed arbitrary point inside a loss 
cavity [15, 16]. The result, not trivially expected, is that in such a condition the 
matter subsystem, because of the cooperation induced by energy loss mechanism, 
may be conditionally guided toward a stationary robust entangled state. The 
renewed interest toward entanglement concept reflects the consolidated belief that 
unfactorizable states of multipartite system provide an unreplaceable applicative 
resource, for example, in the quantum computing research area [25]. However, the 
realization of quantum computation protocols suffers of the difficulty of isolating a 
quantum mechanical system from its environment. In this sense the cited work are 
also aimed at proposing theoretical scheme to bypass decoherence manifestations, 
so taking its place among intense theoretical and experimental research of the last 
few years [15, 16, 26, 27, 28, 29, 30, 31, 32] 2 . 

The paper is structured as follow: in section II I report the principal step 
and approximation leading to the microscopic derivation of the Markovian Master 
Equation, putting in evidence some peculiar properties of it useful in order to 
demonstrate the NuD Theorem N. In section III I solve the markovian master 
equation when T = 0. In section IV I review the applications to old exemplary 
problem and to new, previous unresolved, problem. In section V I try to justify 
the obtained dynamical behaviour in terms of continuous measurement theory. 

2 Citations [15, 16] represent my first results about noise-induced entanglement. Even if they 
are obtained solving the coupled differential equation system of the block vector components 
describing the open system evolution, they can be considered the first NuD theorem application 
because the way in which I obtain the solution is always the same and it find its generalization 
in the NuD theorem. 
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2. Quantum Markovian Master Equation 

It is well know that under the Rotating Wave and the Born-Markov approxima- 
tions the master equation describing the reduced dynamical behavior of a generic 
quantum system linearly coupled to an environment can be put in the form [1] 

Ps (t) = -i[H s + H LS , Ps (t)} + D( Ps (t)), (1) 

where H$ is the hamiltonian describing the free evolution of the isolated system, 

D( P s(t)) =EEvH(^Hps(*)4M 

-\{Ai{u)Ap{u),p s {t)}\ (2) 
Hls = E E S a , p (u)Al(u)A p (u), (3) 

u a,j3 

^M = ^(r ai/3 M-r^H) (4) 

and 

7a,/9M = r Qj/3 (a;) + r^ a (w), (5) 

T a ^{uj) being the one-sided Fourier transforms of the reservoir correlation func- 
tions. Finally we recall that the operators A a (uj) and A^u), we are going to define 
and whose properties we are going to explore, act only in the Hilbert space of the 
system. 

Eq. (1) has been derived under the hypothesis that the interaction hamiltonian 
between the system and the reservoir, in the Schrodinger picture, is given by [1] 

Hi = Y,A a ®B a , (6) 

a 

that is the most general form of the interaction. 

In the above expression A a = (A a )^ and B a = (B a )^ are operators acting 
respectively on the Hilbert space of the system and of the reservoir. The eq. 
(6) can be written in a slightly different form if one decomposes the interaction 
hamiltonian into eigenoperators of the system and reservoir free hamiltonian. 

DEFINITION 1. Supposing the spectrums of Hg and Hb to be discrete (general- 
ization to the continuous case is trivial) let us denote the eigenvalue of Hs (Hb) by 
£ (t]) and the projection operator onto the eigenspace belonging to the eigenvalue 
e (r\) by 11(e) (11(77)). Then we can define the operators: 



A a {u)= E n( £ )^ Q n( £ '), 

£ — 6=U> 



(7) 
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B a (u) = J2 n(r / )B Q n(r ? '). (8) 

From the above definition we immediately deduce the following relations 

[H S , A a {oj)\ = -ujA a (uj), [H B , B a (u)\ = -uB a (u), (9) 

[ff s ,4Ml = +4H and [H B ,Bl{u)\ = +uBl{u). (10) 

An immediate consequence is that the operators A^uj > 0) e A a (uj > 0) raise 
and lower the energy of the system S by the amount Tuo respectively and that the 
corresponding interaction picture operators take the form 

e iHst A a (u)e- iHst = e- iut A a (u), e iHBt B a (u)e~ iHBt = e- iu}t B a (cu), (11) 

e iHst Al(cj)e- iHst = e +iujt Al{u) and e iHBt Bl (u)B~ iHBt = e+^B^ (u). (12) 

Finally we note that 

Al{uj) = A a {-uj) and = B Q (-u). (13) 

Summing eq. (13) over all anergy differences and employing the completeness re- 
lation we get 

£4h = £a*(-u/) = a* and j2 B i^) = J2 B -(- uj ) = B - ( 14 ) 

LU LU LU LU 

The above positions enable us to cast the interaction hamiltonian into the following 
form 

H I= £ A a {u)®B a {u')= £ Ai(u)®Bi(J). (15) 

a,LU,U>' Ct,LU,LU' 

The reason for introducing the eigenoperator decomposition, by virtue of which 
the interaction hamiltonian in the interaction picture can now be written as 

Hj(t)= ]T e-^ + "')U Q (u;)®£ a (u/), (16) 

a,LU,U>' 

is that exploiting the rotating wave approximation, whose microscopic effect is 
to drop the terms for which to ^ —u', is equivalent to the schrodinger picture 
interaction hamiltonian: 

Hj = J2 M") ® B «i-u) = E M") ® 4H- (17) 

a,u> a,ui 
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LEMMA 2. The Rotating Wave Approximation imply the conservation of the free 
energy of the global system, that is 

[H S + H B ,H]=0 (18) 

2.1. Proof 

The necessary condition involved in the previous proposition is equivalent to 
the equation [H$ + Hb, Hi] = we are going to demonstrate. 



[H S + H B ,H] = [H s + Hb,H i } = [H s ,H i } + [H b ,H i } (19) 

= ^[ff s ,4H]04H + ^i Q ( w ) [H B ,Bi(u)} 

= ojA a (uj) <g> B a (-uj) + uA a (uj) (g> B a (-uj) = 0. 

a,ui a,ui 

where we have made use of eq. (9,10) ■ 
LEMMA 3. The detailed balance condition in the thermodynamic limit imply [4] 

la p(uj) = e-^iapi-u) (20) 

where (3 = {ksT)- 1 ■ 

COROLLARY 4. Let us suppose the temperature of the thermal reservoir to be 
the absolute zero, on the ground of Lemma 2 immediately we see that 

7a/3 (u;<0) = ■ (21) 

Let us now cast eq. (1) in a slightly different form splitting the sum over the 
frequency, appearing in eq. (2), in a sum over the positive frequencies and a sum 
over the negative ones so to obtain 



D( PS (t)) 
£ ^{u){A^u)p s A a \uj) 

w>0,o,/3 

- ±{Att(u,)A°(u,), Ps }) 

+ ]T -f^(-u)(A^(u)p s A^u) 

u)>0,a,/3 

- l{A°{ U )Aft(u,),p s }), (22) 
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where we again make use of eq. (13). In the above expression we can recognize 
the first term as responsible of spontaneous and stimulated emission processes, 
while the second one takes into account stimulated absorption, as imposed by the 
lowering and raising properties of A a {uj). Therefore if the reservoir is a thermal 
bath at T = the corollary 4 tell us that the correct dissipator of the Master 
Equation can be obtained by suppressing the stimulated absorption processes in 
eq. (22). 

3. NuD Theorem 

We are now able to solve the markovian master equation when the reservoir 
is in a thermal equilibrium state characterized by T = 0. We will solve a Cauchy 
problem assuming the factorized initial condition to be an eigenoperator of the free 
energy H$ + H B . This hypothesis doesn't condition the generality of the found 
solution being able to extend itself to an arbitrary initial condition because of the 
linearity of the markovian master equation 1 . 

NUD THEOREM 5. // eq. (1) is the markovian master equation describing the 
dynamical evolution of a open quantum system S, coupled to an environment B, 
assumed to be in the detailed-balance thermal equilibrium state characterized by a 
temperature T=0, and if the global system is initially prepared in a state p(Q) = 
PB (0)ps(0) so that {H s + H B )p(0)(H s + H B ) = E 2 L p{0), where E L = E s + E B 
is the free energy of the global system then ps(t) is in the form of a Piecewise 
Deterministic Process [lj, that is a process obtained combining a deterministic 
time- evolution with a jump process. 

3.1. Proof 

The weak-coupling assumption is equivalent to p(t) = p B (0)ps(t). The above 
equation can be used to derive the reduced density matrix p(t) tracing over the 
environment degree of freedom. Let us choose a factorized base B in the tensor 
product Hilbert space made of eigenvectors of Hs and H B 

B = {\E b ,X Eb >\E s ,\e s >}, (23) 

where {E B } and {Es} define respectively the spectra of H B and Hs and {Xe b } 
and {Xe s } their relative degenerations. Let us remember that we have made 
the semplificative hypotheses of discreteness of {E B } and {Es}. In addition we 

It is out of relevance to consider initial condition having non-zero coherence between the 
environment and the system because it is not possible to resolve them in the reduced dynamics 
obtained tracing on the environment degrees of freedom. 
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assume, also for easyness, that {Es} is bounded from below and made of isolated 
points. On the ground of these chioses the total density matrix can be written as 

P{t)= E p(E B ,X EB ,Es,X Es ,E' B ,X' EB ,E' s ,X' Es ,t) (24) 

Eb ,Xe b ,Es,Xe s ,E' b , X' Eb ,E' s ,X' Es 

\Eb,Xe b > \E s ,Xe s >< E' b ,X' Eb \ < E' s ,X' Es \. 

The Lemma 2. imposes a strong selection rule on the indices of the summation, 
that is : 

E s + E B = E' s + E' B (25) 

by virtue of which the trace over the degrees of freedom of the environment, that 
can be written as 

Ps(t)= £ p(Es,XEs,Es,X' Es ,t)\E s ,XE s ><E's,X' Es \, (26) 

E s ,X Es ,E' s ,X' Es 

is immediately obtained: 
Ps(*) = E( £ p(Es,XE s ,X'E s ,t)\Es,XE s ><Es,X' Es \) = J2PEs(t)- (27) 

E s Xe s ,X' Es E s 

The ensemble {Es} can be put in biunivocal correspondence with the natural 
ensemble so that 

N 

Ps(t) = J2Mt), (28) 

i=0 

where N is the natural index corresponding to the maximum of {Es}. The oper- 
ators pi(t) satisfy a big orthogonality condition that is 

Pi (t) ■ Pj (t) = pfSij. (29) 

The last two equations demonstrate that the evolution is a piecewise process (PP) 
or, equivalently, a statistical mixture of alternative generalized trajectories. These 
last are generalized respect to F.Petruccione and H.J.Carmochael approach, which 
leads to ps(t) = J2i IV'i >< V'il- The last expansion, in terms of proper trajectories, 
is obtainable from the mine if and only if we are able to dygonalize the spectral 
correlation tensor, that is known to be always possible, but nobody is able to do 
it, with exception of few highly symmetrical systems. In order to demonstrate the 
formal equivalence between the two approaches we have to demonstrate that the 
alternative processes are deterministic (PDP) or, equivalently, that every of them 
is representable in the form of an evolutionary equation. 



On the ground of Lemma 2.1. and its corollary the markovian master equation 
at T = can be written as 

p s (t) = -i[H s + H LS , Ps (t)} + D(p s (t)), (30) 



D( PS (t)) = e E wH(a»Hp5(*)4h 

-l{4(u,)^H,p5(*)}). (31) 
Let us now substitute eq. (28) into eq. (30), so obtaining 

N ■ N 

= E (32) 

i=0 i=0 i,uj>0,a,f3 

- l -{AP\u)A«(u), Pi }). 

N being the natural index corresponding to the initial eigenvalue of Hg 2 ■ Let us 
observe that 

E 7 a '^M^V)Po^ at M = (33) 

LO>0,a,f3 

so that in the second summation the index i starts from 1 and then the eq. (32) 
can be written as 



ES 1 \fii(t) + -[H ,pi(t)] 

+ E 7 a ' /3 M(^(^)p m ^t( w) + I { A' 3 t( w) ^( W ),p l })] 

ui>0,a,/3 

+ lp- N (t) + UH ,p N (t)}+ E 7 Q '^{^ t H^ a H,ft}]=0. (34) 

12 u»0,a,/3 Z 

Since the big orthogonality of {p{\ the addenda in the above equation act on 
disjointed subspaces of the Hilbert space of H$- This property implies that the 
above equation is verified if and only if the below system of differential coupled 
equations holds 



2 At T=0 only spontaneous emission processes are involved in the dynamics, so that the initial 
eigenvalue of Hs is the maximum permitted eigenvalue. 
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Pi(t) = ~[H , Pi (t)] (35) 
+ l a ' P {u){AP{u)p i+1 A«\u)-UAP\u)A«{u), Pi }). 

uj>0,a,/3 

It is immediate verify by substitution that 

Pl (t) = U(t)f l (t)U\t), (36) 

where, in particular, 

f N (t)= PN (0) (37) 
and U(t) = e~i Bt , U*(t) = e* Bh , B being 

B = Ho ~in S T^M^M^M = Ho - ±H', (38) 

w>0,o,/3 

with H' hermitian. Finally, 

fN-j (t)= £ E - E 7 a/3 (^)7 Q ' /3 '(^)7 a " /3 "K)-7 aJ/3J (^) (39) 

w',a',l3'w",a",l3" ,a' 

x /* C C ... f dt'dt\..dt j U-\t')AP\u')U(t')U-\t")AP" (u n )U(t n )... 
Jo Jo Jo Jo 

U-\t j )A^ j {J )U(t j )f N (t j ) C/t (t j )A aJ t (^' ) C/t- 1 (f ) . . . 

£/t(i" )A Q "t(o;")C/t-i(f' )[/t(t')A a 't(o; , )^ t " 1 (t / ) ) j = l,...,N m 

This concludes the proof and, in addition, ensures that the dynamical pro- 
cesses, whose statistical mixture gives the open system stocastic evolution, are 
deterministic. This demonstrates that the evolution is representable as a Piece- 
wise Deterministic Process (PDP) [1]. The found solution generalizes the PDPs 
introduced by H.J.Carmichael and formalized by F.Petruccione and H.P.Breuer. 
Actually, it is applicable also when the Markovian Master Equation isn't in the 
Lindblad form. This, as already highlighted, in general, introduces simplification 
in the further calculations, but because of the difficulty to recast the equation in 
this form the results obtained are in general merely formal. Tough the eq. (39) 
seems complicated to use it is a powerful predictive tool. I have tested it in a lot of 
contests. Two of my results are already published, others will be object of future 
papers. All of them contain the same (predictable?) result: multipartite systems, 
discarding the physical nature of the parts and of the environment, can exhibit 
entangled stationary states towards the system can be guided by a probabilistic 
scheme of measurement. 



11 



4. Applications 

The best test of a new mathematical tool is its ability to reproduce old re- 
sults and to predict new results. In the contest of old exemplary problem I have 
reproduced: 

- the photocounting formula [33, 8] 

- the enviroment-induced entanglement between two two-level not-direct-interacting 
atoms placed in fixed arbitrary points in the free space [33, 34, 35, 36] 

- Carmichael unravelling of the master equation [33, 5] 1 

Moreover I have tested the NuD theorem's predictive capability solving the 
dynamics of: 

- two two-level dipole-dipole interacting atoms placed in fixed arbitrary points 
inside a single mode cavity in presence of atomic spontaneous emission and cavity 
losses [16]. 

- N two-level not-direct-interacting atoms placed in fixed arbitrary points inside 
a single mode cavity in presence of atomic spontaneous emission and cavity losses 
[15]. 

- A bipartite hybrid model, known as Jaynes-Cummings model, constituted by 
an atom and a single mode cavity linearly coupled and spontaneously emitting in 
the same environment (work in progress) 

- Two harmonic oscillator linearly coupled and spontaneously emitting in the 
same environment (work in progress) 

5. Was the dynamical behaviour of markovian quantum systems 
predictable? The open system in an open-observed environment model 

The physical basis to answer this question is provided by continuous measure- 
ment theory by virtue of which we can regard the Markovian Master Equation 
as the evolutionary equation describing an open quantum system, whose environ- 
ment is continuously monitored. Indirect quantum measurements scheme (Bragin- 
sky and Khalili, 1992, [14]) seems to go to the right direction because it theorizes 
the solution of dynamical problem, here analyzed, as statistical mixture of inde- 
pendent events referred to as non-selective measurement. An indirect quantum 
measurement can be viewed as consisting of three elements. The first element is 
the quantum system of interest and it is called quantum object. The second ele- 
ment is the so-called quantum probe on which the measurement is performed. The 

1 Carmichael interpretation of unravelling is made in a slight different way. He evidences 
the possibility of different unravelling about which I'm trying to better understand. Instead 
F.Petruccione and H.P.Breuer emphasize the non selective-level of measurement that imply an 
independent -from a particular choice of the environment measured observable-, unique unravel- 
ling. This way to think go to my same direction. 
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third element of the scheme is a classical apparatus by which a measurement on 
the quantum probe is performed. 

If we identify the quantum probe as the environment and replace the third 
element with a quantum measurement ideal device (characterized by its own Hilbert 
space) able to detect all the environment elementary excitations coming from the 
system and to make a click as evidence of the detection, it is possible to show that 
the Markovian Master Equation is also obtainable from the microscopic model 
describing this ideal full-quantum measurement scheme. The principal effect of 
this quantum measurement ideal device is to induce loosing of memory in the 
environment respect to the interaction with the system: every time the system 
give an excitation to the environment the device capture it. I name it Memory- 
Cleaner (M-C). 

Let us suppose that a t = the environment is in its ground state (thermal 
state at the temperature T = 0), the M-C is in its ground state 1 and the system is 
in an arbitrary excited state and let us ask: how does they evolve? Two composite 
events may happen: 

- the system does not give excitations to the environment; the M-C doesn't 
capture excitations; 

- the system looses one or more than one excitation in the environment; the 
M-C captures excitations. 

In the first case the M-C will not make a click, the environment, the system 
and the M-C resulting undisturbed by their mutual interactions. In the second 
case the M-C will produce one or more than one click as evidence of the fact that 
the system has exchanged energy with the environment and the environment has 
interacted with the M-C. In this second case the state of the system has changed 
because of the interaction with the environment, the state of the M-C has changed 
because of the interaction with the environment, instead the status of environment 
has been reset, respect to the coupling to the system, to its initial state because 
of the presence of M-C device. The time involved in the reset process can be 
considered of the same order of the coherence time characterizing a markovian 
environment. 

This qualitative analysis shows that the introduction of ideal M-C realizes a 
thermal equilibrium condition maintaining the environment temperature constant 
and equal to zero (deviations from equilibrium can happen in a time scale so short 
respect to the time in which the system change appreciably its state that we can 
assume to look at the system evolution in a coarse-grained time scale assuming 
the environment evolution to happen over time which are not resolved). 

From a microscopic point of view this situation is represented by the hamilto- 

lr The characterization of the spectrum of the M-C is not relevant in this contest. One may 
thinks of a photoemissive source coupled to an environment interacting with one or more than 
one ideal photon detector [1, 5]. 
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nian 

H = H M ~c + H s + H B + Hj + Rbm-c (40) 

containing the free energy of the quantistic M-C (Hm-c)i °f t ne system (Hs) 
and of the environment (H B ) and the interaction hamiltonians (Hi and Hb,m~c) 
describing respectively the coupling between the M-C and the environment and the 
coupling between the system and the environment. The reduced density matrix 
of the system is obtainable by tracing over the environment and M-C degrees of 
freedom the corresponding Liouville equation 

p(t) = - % -[H,p(t)) (41) 

where p(t) is the total density matrix belonging to the Hilbert space given by the 
tensor product of the Hilbert space of the system, of the environment and of the 
M-C. 

The thermal equilibrium condition in the environment, realized by the M-C, 
ensure that 

p(t) = p B (0)ps, M -c(t) (42) 

where /9b (0) is the initial Gibbs state (T = 0) of the environment and ps,M-c is 
the density matrix describing the quantum evolution in the tensor product Hilbert 
space of the system and the M-C. The trace over the M-C degrees of freedom 
{Tr MC {-}) give 

PB(0) Ps (t) = ~[H s ,p B (0)p s (t)] - ^[H B ,p B (0)p s (t)] (43) 
- l -[H IlPB ^)ps{t)] - l jTr M c{[HB,M-CP{t)]} 
and then the trace over environment degrees of freedom (Tr B {-}) give 

ps{t) = ~[H 3 ,ps{t)\ - '-TrB^H^pB^psm (44) 

-^Tr B {Tr M c{[H B , M -c,p(t)]}}. 

The last term in the above equation is zero. This equation has to describes a 
quantum system S coupled to a memory-cleaned environment B. The environment 
characterization ensures that ps(t) satisfies the eq. (30) known as Markovian Mas- 
ter Equation. Actually, this equation takes into account both, the quantum M-C 
device and the environment through the Markov and the Born approximations 
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that can be seen as the memory-cleaned environment approximation and the en- 
vironment equilibrium state assumption. Moreover, the introduction of the ideal 
M-C, beside to give a full justification of Born-Markov approximation, induces to 
think that an almost memory-cleaned environment is modeled by a not-ideal M-C 
which limit case, a full memory environment, is modeled by a not-working M-C, 
a device unable to capture excitations (to interact with environment) or, equiva- 
lently, by a closed microscopic model containing only the system, the environment 
and their interaction. In this way exact open system dynamics can be obtained 
by eq. (40) just removing the M-C and its interaction. On the ground of this 
considerations it is possible to divide the eq. (30) in three terms: 

- we can recognize in the sandwich terms 

]T ^{uj)A^{uj) Pi {t)A a \uj) (45) 

w>0,a,/3 

the jump terms due to the M-C device and describing the capture of an excita- 
tion of frequency to and the consequent projection (jump) of the vector state into 
the subspace characterized by a free energy lower of u> than the previous of the 
measurement act. 

- we can recognize in the free energy constant terms 

~\ E -f^{A^(Lj)A a ( U ),fH(t)}) (46) 

ui>0,a,/3 

the effective dissipative part responsible of going to zero of populations and co- 
herences of the excited states of the open system of interest and due to the open 
observed environment we are considering. In fact these terms introduce an imag- 
inary frequency in the free-energy spectrum of Hs, so generating an exponential 
decaying of excited free eigenstates. 

- we can recognize in the free energy Lamb-Shift terms 

E S°>P[Aft(u)A a (u), Pi (t)] (47) 

ui,a,/3 

the environment-induced multipartite cooperation not- vanished by the loss of mem- 
ory of the environment. 

The found solution (NuD theorem) tell us that the state of the system is a 
statistical mixture of the free energy system eigenoperators. This fact depends and 
it is consistent with the existence of the M-C measurement device because the act of 
measurement introduces a stochastic variable respect to which we can only predict 
the probability to have one or another of the possible alternative measures. These 
probabilities can be regarded as the weight of the possible alternative generalized 
trajectories and, analytically, they are given by the partial traces of the pi. With 
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this approach the dynamics has to be depicted as a statistical mixture of this 
alternative generalized trajectories. 

Moreover the found trajectories evolve in time in a deterministic way: for 
example the trajectory relative to the initially excited system state is a shifted 
free evolution characterized by complex frequencies that means an exponential 
decaying free evolution. This statement may give the sensation that every system 
has to decay in its ground state because of the observed dynamics. It is in general 
not true. Actually, if the system is multipartite, it is possible that it admits excited 
and entangled equilibrium Decoherence Free Subspace (DFS) (so as it happens in a 
lot of known models), constituted by states on which the action of Hi is identically 
zero and then, if the system, during evolution, passes through one of these states, 
the successive dynamics will be decoupled from the environment evolution. An 
equilibrium condition is reached in which entanglement is embedded in the system 
[15, 16]. 

What could ensure, for example, that an entangled decoherence-free state, if 
existing, has been generated? (For example) the number of click we hear in a 
period of time long enough respect spontaneous emission rate. Actually, if the 
numbers of the clicks is less than the numbers of initial excitations then we can 
say that it has been generated a decoherence-free state [15, 16, 34]. 

If, on the contrary, the system is monopartite it is possible to demonstrate that 
the only possible DFS is generated by the ground state of the system so that a 
monopartite system will loose its internal coherence and the population of excited 
states because of the measurement and the consequent interaction with a memory- 
cleaned environment, unable to induce cooperation among the parts: there exists 
only one part. 

In this case the numbers of clicks has to be the same of initial system excita- 
tions. The M-C sounds as a quantum logic counter. 
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